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The possible transport of fibers by fluid flow in fractures is investigated experimentally in trans- 
parent models using flexible polyester thread (mean diameter 280/im) and Newtonian and shear 
thinning fluids. In the case of smooth parallel walls, flbers of finite length £ — 20 — 150 mm move at 
a constant velocity of the order of the maximum fluid velocity in the aperture. In contrast, for flbers 
lying initially at the inlet side of the model and dragged by the flow inside it, the velocity increases 
with the depth of penetration (this results from the lower velocity - and drag - in the inlet part). 
In both cases, the friction of the fiber with the smooth walls is weak. For rough self-affine walls 
and a continuous gradient of the local mean aperture transverse to the flow, transport of the fibers 
by a water flow is only possible in the region of larger aperture (a > 1.1mm) and is of "stop and 
go" type at low velocities. Time dependent distorsions of the fiber are also often observed. When 
water is replaced by a shear thinning polymer solution, the flbers move faster and continuously in 
high aperture regions and their friction with the walls is reduced. Fiber transport becomes also 
possible in narrower regions where irreversible pinning occurred for water. In a third rough model 
with no global aperture gradient but with rough walls and a channelization parallel to the mean 
flow, flber transport was only possible in shear-thinning flows and pinning and entanglement effects 
were studied. 

PACS numbers: 



I. INTRODUCTION 

Fiber transport by flowing fluids is of interest in many 
areas of physics, biology and engineering: examples in- 
clude the transport of paper pulp [l|[, the manufac- 
turing of fiber-reinforced composites the rheology 
of biological polymers [3] and the motility of micro- 
organisms ^, i5j]. Recently, also, using long optical fibers 
has been suggested as a method to realize distributed in- 
situ measurements (temperature for instance) on natural 
water flows @: this raises the problem of the possible 
transport of the fiber by the moving fluid which often 
occurs in geometries confined by solid walls. More pre- 
cisely, little work has been devoted to flow channels with 
rough walls such as fractures of natural rocks and of ma- 
terials of industrial interest in civil, environmental and 
petroleum engineering. In that case, the interactions of 
the fibers with the walls are particularly important and 
may lead to blockage of the motion of the fibers and/or 
to clogging of the channels. In addition to their practical 
applications, these processes raise fundamental questions 
regarding the motion of flexible solid bodies in complex 
flow fields. 

The objective of the present experimental study is to 
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understand the transport of fibers by the flow and to 
investigate the role of the flow geometry and fluid rhe- 
ology. Here, single fractures are modeled by the space 
(saturated by a flowing fluid) between either two parallel 
plane walls or two complementary rough self-affine walls 
with a relative shear displacement from their contact po- 
sition. This latter configuration allows one to reproduce 
preferential flow channels 0, [1] which are a widespread 
feature of natural fractures and influence strongly their 
transport properties. The walls of the fracture are trans- 
parent to allow for optical observations of the motion of 
the fibers. 

In addition to hydrodynamic forces on the fibers due 
to the relative velocity with the fluid, their motion and 
deformation is influenced by different effects. A first one 
is the interaction forces with the walls: they are par- 
ticularly important when the diameter of the fibers is 
comparable to the local channel aperture or when they 
are close to one of the walls ^ JOi]. Second, tension 
forces reflecting the mechanical cohesion of the fiber are 
present all along its length so that the motion of each 
region influences the other ones. Tension and hydrody- 
namic forces add-up and their spatial variations in shear 
flows or flows with curved streamlines deform the flbers 
(although their length remains constant): this may fi- 
nally lead to entanglement and blockage. Next, these 
deformations are opposed by elastic forces reflecting the 
non-zero stiffness of the flbers: their relative magnitude 
with respect to the hydrodynamic forces is a particularly 
important parameter of the problem [ll[. On the one 
hand, very flexible fibers would seem to be able to follow 
the streamlines but, on the other hand, loops may appear 
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easily and lead to trapping in narrow zones . Another 
key element is the rheology of the fluids which influences 
strongly the hydrodynamic forces on slender bodies, even 
in simple shear flows Finally, inertia (particularly 

in regions of large spatial flow velocity variations) and 
gravity may also be of importance by inducing motions 
transverse to the flow and towards the walls. 

In the following, the feasibility of fiber transport and 
its dependence on the mechanisms discussed above are 
studied in three different model fractures. A first one 
has smooth walls and is used as a reference case and the 
two others have rough walls with self affine geometries. 
For one of the rough fractures {F3) , the mean planes of 
the walls are parallel. For the other one (-^2), they have 
a small angle resulting in a non-zero transverse gradient 
of the mean aperture: this wedge shape mimicks the edge 
of many natural fractures. In this work, the infiuence of 
the geometry of the fibers on their transport is analyzed 
by comparing the motion of finite length segments (still 
with an aspect ratio greater than 200) and of continu- 
ous threads injected at the inlet of the fracture. The 
influence of deformations due to flow velocity gradients 
could be investigated by using flexible thin flbers made 
of polyester thread. Special attention has been brought 
to the mechanisms of pinning during the motion of the 
fibers and of possible depinning after enough elastic en- 
ergy has been accumulated (see sections IlII Bl and IIV Ap . 
Finally, the influence of the fluid rheology has been stud- 
ied by comparing fiber transport by Newtonian and shear 
thinning fluids. 



II. EXPERIMENTAL SET-UP AND 
PROCEDURE 

A. Fiber characteristics 

The experimental fibers are prepared from commer- 
cial polyester thread used for needlework and made of 
two strands twisted together. The section is not circular 
so that the diameter varies between 220 and 340 /xm. 
The specific mass of the fiber is p — 1.8±0.1 x 
10~^ kg/ . Its bending stiffness EI (ratio of the ap- 
plied bending momentum by the curvature) is of the or- 
der of 10~^ kg.m^/s'^ (a value similar to that reported 
in ref. [ij] for a comparable material). This value has 
been estimated by measuring the deflection under its 
own weight of a horizontal fiber segment attached at one 
end [13. 

The choice of these multifilament fibers refiects a trade- 
off between several requirements. The fiber diameter is 
large enough to be clearly visible and its flexibility is 
both high enough to allow for deformations by the fluid 
velocity gradients and low enough so that coils do not 
build up too easily (see introduction). We verified indeed 
that monofilament wires of similar (and even smaller) 
diameters were too rigid and retain often, in addition, a 
permanent curvature. 



Both segments of fiber and continuous flbers were 
used in the present work. The segments are cut out of 
polyester thread with a length 20 mm < £ < 150 mm. 
The continuous flbers are also cut out of a spooled thread 
but with a length larger than that of the fracture: they 
are progressively injected from above into the fracture 
under zero applied tension conditions. A speciflc proce- 
dure is used to avoid blunting the ends of the fibers and 
the flbers are carefully saturated with liquid prior to the 
experiments. 



B. Model fractures 



I 
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FIG. 1: Scliematic view of tlie experimental models, (a) frac- 
ture with flat parallel walls - (b) fracture with complementary 
self-afiine walls with a relative displacement u. For all cells 
: w = 90 mm, at — 20 mm, = 5 mm, = 52 mm, 
h ~ 20 mm and L — 288 mm. The mean flow is vertical and 
paraUel to x 

The fracture models are manufactured with the same 
technique as in Ref. [16^ by carving two parallelepipedic 
plexiglas blocks using a computer controlled milling ma- 
chine. The two blocks are then clamped together in a 
preset position determined by the geometry of the sides 
of the block: these act as spacers leaving a controlled 
interval between the surfaces for the fluid flow (Fig. [1]). 
This procedure allows for the realization of model frac- 
tures with rough (or smooth) walls of arbitrary geome- 
tries and three of them have been used. 

The first sample (refered to as Fl) has smooth par- 
allel plane walls separated by a fixed distance a{x, y) = 
0.65 mm (Fig. [T^). The two other models have comple- 
mentary rough walls with a self affine geometry of char- 
acteristic Hurst exponent H = 0.8: this value is similar 
to that measured on granite fractures (l7j and in vari- 
ous other materials [18]. The amplitude of the rough- 
ness, represented by the difference between the heights 
at the highest and lowest points of the surface is equal 
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to 26 mm: this corresponds also to values measured on 
natural rock fractures |19|]. For these two models, the 
wall surfaces are perfectly complementary: in order to 
generate aperture fluctuations, they are positioned with 
a small relative shear displacement u — 0.65 mm in the 
direction y normal to the flow (Fig. [Da). 

The resulting aperture maps are displayed for the two 
models in Figures HI and [TUl using a color code. The bot- 
tom curves represents the transverse profile in the y di- 
rection of the average along x of the aperture. In model 
there is a small angle between the mean plane of 
the surfaces mimicking the edge of natural fractures. As 
a result, the average a of the aperture along the direc- 
tion X of the flow dec reases wit h the coordinate y from 
a{y = 0) = 1.1 mm to a{y — w) — 0.8 mm with a global 
mean value a = 0.91 mm (Fig. Uj;); w is the width of 
the cell and y = and y = w correspond to the sides of 
the experimental region. Note that the wedge shape of 
the model results in a global color gradient in the direc- 
tion y in the maps of Fi g. H^ -b. The root mean square 
amplitude aa{y) = (< {a{y) - a{x,y)Y >x)^^'^ of the 
fluctuations of the apert ure i n the direction x with re- 
spect to the mean value a{y) =< a{x, y) >x is equal to 
^a{y) = 0.122 mm and independent of y. Therefore, the 
fluctuations remain of constant absolute amplitude while 
the mean aperture varies with y. 

In model FS, the mean planes of the fracture walls are 
parallel with a distance a ~ 0.65 mm comparable to the 
value for model Fl. The local fluctuations of the aperture 
are of the same order of magnitude as the deviations from 
the linear trend in model F2 (Fig. [TUh). No global trend 
in the color shade at the scale of the global size of the 
fracture is visible in the map of Fig. [TOk-b. Comparing 
experiments realized with the two systems allows one to 
discriminate between the influence of the roughness and 
that of the global transverse aperture gradient. 



C. Experimental set-up and procedure 

The model fractures are held vertically; liquid is sucked 
uniformly at the bottom side at a constant flow rate and 
reinjected at the top into an open bath of area 5 x 100 mm 
and depth l}, ~ 10 mm covering the inlet of the model. 
This design allows one to introduce the fibers in the bath 
through its open surface and then into the fracture. In 
order to make the injection of the fibers easier, the up- 
per section has a funnel-like "Y" shape, i.e. the mean 
aperture of the fracture increases with height in the top 
52 mm of the model (distance le in Figs.flji-b). 

The model is illuminated from behind by a light panel 
and a digital camera provides images with 1024 x 768 
pixels at a rate of 30 frames per second and with an 
exposure time of 1/300 s. The length of the field of view 
is 150 mm in the vertical direction parallel to the flow and 
the top of the aperture maps corresponds to a distance 
of 110 mm from the inlet of the fracture. 

On each picture, the location and geometry of the fiber 



is determined by a binary thresholding technique. For 
most experiments, except in some cases in which the fiber 
gets pinned, its distortions are moderate and it remains 
overall aligned with the direction of the flow. Therefore, 
at each distance x from the inlet, there is only one point 
along the fiber centerline, and, therefore, only one value 
of y. ^From these data, spatiotemporal-diagrams like 
those displayed in Figs. [3k-b are obtained. The vertical 
scale corresponds to the vertical distance x, the hori- 
zontal one to time and the color code corresponds to the 
distance y at the corresponding value of x. Qualitatively, 
these diagrams provide information on the global motion 
of the fibers, on their rotation and on their deformation. 
For instance, diagrams with uniform colors on horizon- 
tal lines correspond to a pure translation of the fiber in 
the direction parallel to the mean flow; in contrast, color 
variations of these lines imply lateral motions. Quanti- 
tatively, the velocity of segments of fibers with a finite 
length is take equal the variation between two successive 
images of the vertical distance Xm of their center of mass 
to the inlet side of the model. For continuous fibers, one 
uses instead the variation of the distance Xt of the tip 
of the fiber to the inlet. In both cases, the vertical ve- 
locity component of the center of mass (respectively tip) 
of the fiber will be referred to as V/. Continuous fibers 
are more easily retrieved after the experiments than fiber 
segments when trapping occurs: for that reason, most 
experiments in rough fractures have been realized with 
continuous fibers in order to obtain more data. 



D. Characteristics of the fluids 

The Newtonian fluid used in the experiments is high 
purity water (Millipore - Milli-Q grade) with a specific 
mass p = 10"^ kg/m'^ and a dynamic viscosity ^ao — 
1 mPa.s. The shear thinning fluid is a solution of high 
molecular weight Scleroglucan (Sanofi Bioindustries) in 
water at a concentration Cp = 1000 ppni. The variation 
of the effective viscosity ji with the shear rate 7 (see 
Figure 2 of ref. [l^]) is well described by the Carreau 
function: 



(1 + (7/70)2)- 



-(MO - Moo) + Mo 



(1) 



with n = 0.26, 70 = 0.026 ± 0.004 s-\ no = 4500 ± 
340 mPa.s. A shear thinning domain of shear rates in 
which p oc 7("~i) is bounded by two Newtonian domains. 
At low shear rates (7 <C 7'o), the viscosity is constant 
with /i = /iQ. At very high shear rates (7 ^ 700 = 
7o(mo/Moo)^^^^"), the effective viscosity should reach the 
limiting value ^00 ■ Practically, and following theoretical 
expectations, is taken equal to the viscosity of the 
solvent {i.e. fioo = 1 mPa.s for water) because it cannot 
be measured directly (the corresponding transition shear 
rate 700 — 2300 s~^ is above the useful range of our 
rheological measurements (7 < 100 s~^)). 
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E. Characteristic velocities and shear rates in the 
experiments 

In the present experiments, the mean flow velocity 
ranges from U — 50 to 400 mm/s. For water, the corre- 
sponding Reynolds number defined as Re = Uap/ varies 
from 40 and 320: substantial inertial effects are therefore 
expected (21j . and will be larger for rough fractures. An 
alternative value of the Reynolds number is obtained by 
using the diameter of the fiber and its relative velocity 
with respect to the fluid; it remains however higher than 
1, confirming the need to take into account inertial ef- 
fects. 

Regarding the shear rate 7, it is equal to zero in the 
center part of the aperture and reaches a maximum at 
the walls. For a Newtonian fluid, 7^, — 6U/a (neglecting 
the influence of the fiber and of the roughness of the 
walls). For shear thinning fluids [13], jw depends on the 
relative values of the shear-rate 7 and of the parameters 
70 and 7oo (see Sec. Ill D|) in the flow. At very low mean 
flow velocities U for which 7^, < 70, the viscosity would 
be equal to /iq in the full fluid volume and the results 
would be the same as for a Newtonian fluid. Using the 
numerical values from Sec. IIIDl this would require that 
U < a7o/6 = 2.16 x 10"'^ mm/s, a value well below the 
present experimental range [U > 50 mm/s). The fluid 
has therefore a shear thinning behavior in a part of the 
flow volume. 

In the opposite high velocity limit, a second Newtonian 
layer with p = develops at the walls if 7^, > 700. As- 
suming a truncated power law rheological curve (see [23|). 
the corresponding threshold value of U may be estimated 



as: U 



l/n 



3n 



2n+l 



Using the numerical values 

from Sec. IIIDi this leads to U — 2 mm/s, again below 
the experimental range of velocities used in the present 
work. A second Newtonian layer of viscosity /Jac appears 
therefore on the walls in our experiments and the effec- 
tive local viscosity increases from at the walls to ^0 
in the center of the gap. The thickness Sz of the New- 
tonian layer of viscosity fico on the walls at the highest 
experimental ffow velocity U = 400 mm/s may be esti- 
mated from the relation 2Sz/a = 1 — ajoo/ (pU) valid if 
iw ^ 7oo : the combined thicknesses of these two Newto- 
nian layers on both walls are of the order of 52% of the 
gap a. 



III. FIBER TRANSPORT BY WATER IN 
SMOOTH AND ROUGH MODEL FRACTURES 

A. Fracture with smooth walls (model Fl) 

1. Transport of fiber segments 

As a first step, reference experiments were carried out 
on the model fracture Fl with flat smooth walls, starting 
with fiber segments of length ^ ranging from 20 to 150 mm 
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FIG. 2: Variation of the normalized velocities Vf/U of the 
tip of continuous fibers as a function of the distance xt of the 
tip from the injection side of the model: (■) U — 100 mm/s, 
(•) U = 150 mm/s, (a) U = 290 mm/s. Variation of the 
velocity dxm/dt of the center of mass of fiber segments of 
length I — 20mm as a function of Xm'- (□) U = 100 mm/s, (o) 
U = 150 mm/s, (A) U = 290 mm/s. Dashed line : Vf/U = 
1.1 — 105/a:t. Inset: variation of Vf/U as a function of 1/xt 
for the continuous fiber (same symbols as in main graph). 
Dash-dotted Une: estimated value Vf/U — 1.35 from ref. ^ 



and at mean velocities U of the fluid across the aperture: 
50 < U < 400 mm/s. It should be noted, in addition, 
that few experimental and even numerical results (out- 
side reference Q) are available on the transport of fibers, 
even in the simple geometry of parallel plane walls). The 
present study is centered on the transport of fibers with 
their mean axis parallel on the average to the direction of 
the fiow within ±10°. For each set of values of i and U, 
more than 20 transport experiments were realized. In all 
cases, no deformation of the fibers is visible during the 
motion and they reach a velocity which is constant with 
the distance Xm of their center of mass to the injection 
side within ±10%. This latter result is visible in Figure[2] 
in which the ratio of the velocity of the fiber (V/) and of 
the mean fluid velocity (U) is plotted as a function of Xm 
for 3 different U values (open symbols) in a range of 1 
to 3. Taking into account the intrinsic dispersion due to 
the injection process and the variability of the shape of 
the fiber, the variations of Vf/U (< 25%) are too small 
to suggest any definite dependence otVf/U on U. More 
generally, in other experiments performed at different ve- 
locities and for fibers of different lengths, Vf/U always 
ranges between 1 and 1.35 at all values of Xm- 

These observations can be compared to the numeri- 
cal simulations by Sugihara-Seki 9] of the longitudinal 
velocity of a neutrally buoyant circular cylinder located 
midway between the walls of a plane channel. For a ratio 
of the diameter of the cylinder to the distance between 
the walls equal to 0.5 as in the present experiments, the 
normalized predicted velocity Vc/U is equal to 1.35 0]; 
this is consistent with the upper bound of the present ex- 
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perimental results (dash-dotted line in Fig. [2]). The lower 
values would then correspond to cases in which the fiber 
is not in the center of the interval or in which it has a 
wavy shape which obstructs partly the flow field. 




FIG. 3: Spatiotemporal diagram of the transport of a contin- 
uous fiber by a downward flow of water in two model frac- 
tures. Horizontal axis: time - vertical axis: distances in the 
vertical direction counted from the top of the model (see 
Fig. [T|. Color code: transverse coordinate y{x) of fiber at 
the corresponding distance x. (a) Fracture Fl with smooth 
walls. {U = 120mm/s). y ranges from 52.9mm (green shade) 
to 63.3 mm (light blue) - (b) Fracture F2 with rough walls 
{U = 160mm/s). y ranges from 58.4mm (green) to 71.5mm 
(dark blue). End of fiber inserted initially in the high perme- 
ability region of the fracture (a > 1 mm). 



2. Transport of continuous fibers 

The above results are now compared to the case of a 
continuous fiber lying initially at the top of the model 
with no external force applied and dragged by the flow- 
ing fluid into the space between the fracture walls once 
its tip has been inserted into it. Figure [3^ displays a spa- 
tiotemporal diagram obtained for such an experiment by 
the technique described in Sec. Ill CI At a given time t, 
the color is constant along x in the colored zone, indicat- 
ing that the fiber remains straight and vertical. The color 
varies slightly with time which reflects a small global 
transverse displacement. The boundary of the colored 
zone marks the coordinate Xt of the tip at the corre- 
sponding time t: its slope (proportional to the velocity 
Vf) increases with t. Quantitatively, Fig. [5] shows (dark 
symbols) that the normalized velocity Vf /U of the tip of 
the fiber increases with Xt (U is the mean flow velocity 
in the constant aperture part of the model) . 

We show now by a simple model that this variation of 
Vf/U reflects the lower flow velocity in the upper 52 mm 
of the fracture due to the larger cross section of the flow 
(this upper part has a "1^" shape with an aperture in- 
creasing towards the top). We assume that the flber: a) 



remains straight and vertical so that the velocity of all 
its points is equal to Vf and: b) is located in the center of 
the gap where the velocity has a maximum value Vm{x) 
(assumed to be independent of the transverse coordinate 

y)- 

As a very first approximation, the local vertical viscous 
force fv per unit length on the fiber is taken equal to /„ = 
/iC(VM(a;) — Vf) in which C is a geometrical constant. 
Due to the low angle and the small distance between 
the walls for x > 0, the flow field is quasi parallel with a 
Poiseuille-like profile between the walls at a given vertical 
coordinate a: > (for a Newtonian fluid and neglecting 
flow perturbation by the fiber). In the upper fluid bath 
{x < 0), this result is not satisfied but the fluid velocity 
is low enough so that its contribution is negligible. 

In order to comply with volume flow rate conservation, 
the mean fiuid velocity V{x) at the distance x where the 
local aperture is a(x) must satisfy: V{x) a{x) — Ua — 
cst{x) in which a is the the aperture in the parallel part 
of the cell, where the mean velocity is U . The total 
viscous force on the continuous fiber between its lower 
tip at the vertical coordinate Xt and the surface of the 
upper bath is then : 



F^{xt) = / U{x)dx = tiC BUa I - [h + Xt)Vf 

a(x) 

(2) 

a{x) — a for x > le and a(x) decreases linearly from to 
a with X for < x < le- As mentioned above, the velocity 
of the fluid in the upper bath (x < 0) has been neglected. 
The coefficient B would be equal to 3/2 for a Poiseuille 
proflle with no obstruction effect. The integral may be 
readily computed and Vf may be estimated by taking 
Fy = 0: this is equivalent to neglecting inertial forces as 
well as friction on the walls which has been found to be 
approximately valid in the case of the previous segments 
(the flber moved roughly at the maximum velocity in the 
gap). This leads to: 



BU 



Xt — le -\ In — 

Oe — a a 



The ratio Vf/U may then be rewritten as: 

u 



Xt + lb 



(3) 



(4) 



in which b* and I* are two constants and Xt + lb is the 
vertical distance of the tip of the flber to the free surface 
of the bath. Using the parameter values given in Fig. [T] 
leads to the theoretical prediction: b* = B = 3/2 and 
I* = lb + le{l — a£n{ae/a)/{ae — a)) ~ 57mm. Exper- 
imentally, b* and I* can be determined conveniently by 
plotting Vf /U as function of l/{x + lb) (inset of Fig. 
The variation is linear as expected and a linear regres- 
sion over the data gives the values b* = 1.1 ± 0.1 and 
I* — 105 ± 15 mm. These values are valid for all three 
velocities U investigated and also lead, as might be ex- 
pected, to a good flt (dotted line) in the main graph of 
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y (mm) 

FIG. 4: (a-b) Snapshots of continuous fibers (in black) dur- 
ing experiments realized in model F2 with two transverse in- 
jection distances (vertical arrows). The images of the fiber 
are superimposed onto a color coded map of the aperture 
field (red: a < 0.4 mm, blue: a > f.4mm). Field of view: 
90 X f 50 mm. Distances in the vertical x direction are counted 
from the top of the model (see Fig. [T]). Vertical arrows mark 
the location of the entry point of the fibers inside the field of 
view, a) [/ = f60mm/s; horizontal arrows: 5 pinning posi- 
tions h) U — 250 mm/s. The fiber is pinned outside the image 
at a point of coordinates {x = f50mm and y — 60mm). (c) 
variation of the aperture averaged along the distance x par- 
allel to the flow as a function of the transverse distance y. 

of Fig. [21 The difference between the experimental and 
theoretical coefficients likely reflects perturbations of the 
flow profile by the fiber so that both B and C may differ 
from the theoretical value and depend on x. Finally, ne- 
glecting the influence of the inlet would only be possible 
for X >> I*. This is not the case here since the global 
fracture length is only ~ 3r. 



B. Fiber transport in a wedged fracture with 
rough walls 

While, for fractures with smooth plane walls, the fluid 
velocity is constant with x and y and no deformation of 
the fibers is observed, the flow fleld is very complex in 
rough fractures with large shear stresses in the fracture 
plane. If the fiber is sufficiently flexible, it may get bent 
by these stresses. Also, friction of the fibers with the 



walls is larger than for the smooth fracture of Sec. IIII Al 
and it screens the influence of the inlet much more effi- 
ciently. All the experiments with rough walls described 
in the rest of the paper have been realized with a contin- 
uous fiber. 

Fiber transport by a water flow has been studied in 
fracture F2 which has a self-affine roughness geometry 
with a small angle between the mean planes of the walls. 
In this model, one has both random variations of the 
velocity due to the roughness and a mean transverse gra- 
dient. Since the mean aperture a(y) varies continuously 
with y (Fig. (31;), this allows us to test the minimum value 
of a{y) (averaged over the length L) for which fiber trans- 
port all along the fracture is possible. In the experiments 
using water as the carrier fluid, this minimum value was 
found to be of the order of 1 mm (typical mean aperture 
in the left side of the maps of Figs. [4]). A view of a fiber 
obtained during such an experiment is superimposed onto 
the aperture field in Figure [3^. In the low aperture part, 
most fibers get pinned in the first cm past the initial 
large aperture section at the top of the model. How- 
ever, in a few experiments performed at high flow rates 
{U > 250 mm/s), the tip of the fiber gets first pinned but 
a loop builds up downfiow of the tip (Fig. [4)d); then, the 
bottom of the loop moves sideways into the high aper- 
ture region where the velocity is largest. Finally, the loop 
pursues its downward motion until it reaches the bottom 
of the model. 

Let us now concentrate on the dynamics of a fiber in- 
jected in the high permeability (large aperture) region at 
the left. Figure [3)3 displays a spatiotemporal diagram of 
such an experiment: it differs strongly from the diagram 
of Fig. [3^ obtained for a model with plane walls. The 
motion of the tip is not continuous but corresponds to a 
sequence of "stop" and "go" phases marked respectively 
by horizontal and oblique sections of the boundary of the 
colored zone in Fig. The duration of individual pin- 
ning events may range from less than 1 s to nearly 10 s 
(about 10 % of the total transit time of the fiber along the 
fracture). One also observes that pinning is always initi- 
ated at the tip at the fiber and not in the zones located 
behind it. 

The sites where the fiber gets pinned are shown by 
horizontal arrows in Fig. 2^. Several (but not all) of them 
are located in regions where the aperture is below 1 mm 
(yellow-green shades in Fig. 0^). This is still twice the 
diameter of the fiber so that pinning cannot result solely 
from blockage in a constriction too small to accomodate 
the fiber. Moreover, the fiber is observed to cross other 
regions of the fracture with a similar aperture without 
getting stopped. Blockage of the fibers may also result 
from a large local slope of the surfaces, even if the local 
spacing remains large. In this case, the hydrodynamic 
forces may be too small to induce enough bending of the 
fibers so that they follow the slope of the surface and do 
not touch it. Examining the geometry of the individual 
wall surfaces around the pinning points tends to confirm 
this hypothesis. Finally, the inertia of the fiber may also 
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favour contacts with the walls and blockage: this effect 
should increase with the fluid velocity. 



100 





10 20 3010 20 30 
y (mm) y (mm) 



10 20 30 
y (mm) 



10 20 30 
y (mm) 



FIG. 5: Variation of the fiber geometry during the fifth pin- 
ning event in Fig. (lowest horizontal arrow). The left- 
most (respectively rightmost) picture displays the shape of 
the fiber immediately after pinning (respectively before de- 
pinning). Time lapse between pictures: 0.f66 s; total dura- 
tion of trapping event: f s. Arrows: location of points of 
maximum curvature referred to as 0, f and 2. 

More information on the pinning-depinning process is 
obtained from the geometry of the fiber. The variations 
of the code color with time and distance in Fig. in- 
dicate that, in contrast to the case of a smooth fracture 
(Fig. [3^), the fiber does not remain straight: it displays a 
buckling instability and meanders of shape and location 
varying with time appear. These variations are visual- 
ized directly in Figure displaying the deformations of 
the fiber. During the pinning event, the tip of the fiber 
remains motionless while three bumps appear behind it; 
their amplitude increases with time while they propagate 
toward the tip. The deformation of the rear part of the 
fiber is much weaker while, in some cases, it gets trans- 
lated sideways. 

We characterized quantitatively the deformation of the 
fibers by the maximum absolute value \R\ of the radius 
of curvature in the three bumps and by the distance Sx 
of the corresponding points to the tip of the fiber (mea- 
sured along axis x). The variation of |i?| and 6x with 
time for all three regions in this same experiment is dis- 
played in Fig. [SI The distance Sx decreases with time 
with a relative velocity of the order of lOmm/s lower 
than the fluid velocity (160mm/s in Fig. 3^) and similar 
for the three points of interest while the corresponding 
radii \R\ decrease strongly. The deformation is more im- 
portant in the two first bumps for which the radius of 
curvature is of the order of 9.5 mm at the time when 
depinning occurs. These deformations reflect the accu- 
mulation of bending elastic energy behind the tip which. 
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FIG. 6: Variation with time of the absolute values |-R| of the 
radius of curvature of the fiber at the 3 maximum curvature 
points designated by arrows in Fig. [5] Inset: time variation 
of the vertical distances Sx between these 3 points and the 
fiber tip. Symbols for point 0: (o,»); for point I: (□,■) and 
for point 2: (A, A). 
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FIG. 7: Variation of the ratio Vf/U of the fiber and fluid ve- 
locity as a function of the distance xt from the inlet of model 
F2 for three experiments with continuous fibers injected at a 
same location, (o) and (•): water fiowing at respective veloc- 
ities U = 160 mm/s and 210 mm/s. (*): fiow of a polymer 
solution at f/ = 220 mm/s. Dotted and solid lines: average 
of fiber velocity over field of view respectively for U = 160 
and 210 mm/s. 



finally, gets large enough to overcome the friction forces 
with the wall. Similar scenarios were observed in all cases 
in which pinning takes place even though the time scales 
varies. 

The pinning/depinning process should depend signifi- 
cantly on the hydrodynamic forces and, therefore, on the 
flow velocity U. This influence is analyzed in Figure [7| 
displaying the variation with distance of the normalized 
velocity Vf/U of the fiber tip for two different values of 
U. At the lower velocity U = 160 mm/s, the fiber dis- 
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plays a stop and go motion already visible qualitatively 
in Fig. [3] corresponding to the same experiment. The 
velocity cancels out at each pinning point and increases 
sharply after the release of the fiber to a value much 
higher than the time average (horizontal dashed line). 
For U — 210mm/s, the variations with time of the fiber 
velocity with respect to the mean value are of smaller am- 
plitude (-100%,+50%) instead of (-100%,+150%). Also, 
the fiber only stops completely once while Vf decreases 
by less than 40% on other pinning sites: these are located 
at the same distances x as for U = 160mm/s. These re- 
sults imply that the influence on pinning of inertia forces 
dragging the fibers towards the walls is small. The nor- 
malized averages Vf/U of the fiber velocity with time 
are equal to 0.2 and 0.28 respectively for U = 160 and 
210mm/s. Even if the influence of the low velocity near 
the pinning sites is removed from this average, Vf/U is 
still only of the order of 0.3 for both U values. This is 
lower by a factor of 3 than the values measured in model 
fractures with smooth plane walls (see Fig. [5]). 



IV. FIBER TRANSPORT BY SHEAR 
THINNING FLUIDS IN ROUGH MODEL 
FRACTURES 

A. Fiber transport in fracture F2 

The influence of rheology has been studied specifically 
by replacing water by a shear thinning solution of char- 
acteristics discussed in section IIIDI in the same model 
F2 as above. In Fig. [71 data obtained with this fluid 
are compared to the variations for water discussed in the 
previous section (the points of injection of the fibers are 
the same and correspond to the high aperture path) . At 
a similar velocity (~ 220mm/s), the ratio Vf/U is larger 
for the shear thinning solution than for water, except at 
short distances. The amplitude of the velocity fluctua- 
tions is also smaller but they still take place at the same 
distances x as the pinning sites observed in low veloc- 
ity water experiments. The motion of the fiber remains 
therefore slightly influenced by the local asperities en- 
countered along its path, although much less than for 
water. 

Globally, using a shear thinning fluid instead of water 
reduces significantly the influence of the roughness of the 
walls on the motion of the fibers: their dynamics may be 
then much more similar to that reported in Sec. lIII Al for 
the plane fracture Fl. In particular, the increase of Vf /U 
with the distance x indicates that the reduced drag forces 
on the parts of the fiber located at the inlet influence the 
global motion of the fiber, like for fracture Fl but unlike 
the experiments of Sec. IIII 5] with water. 

The above results correspond to fibers moving, like in 
Sec. IIII Bl in the high aperture part of the width of the 
wedged fracture F2. However, fibers may propagate even 
in the low aperture part, if the fluid is shear thinning 
while they get blocked or build up side loops with water 
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FIG. 8: Variation of the fiber geometry during a pinning event 
for a polymer flow {U = 230mm/s) in the low aperture part 
of fracture F2. The leftmost (respectively rightmost) picture 
displays the shape of the fiber 0.1 s after pinning (respectively 
just before depinning). The definitions of points 0, 1 and 2 
are the same as in Fig. [5] 

(see Fig. 2]). The corresponding motion is then not con- 
tinuous but of a "stop" and "go" type: the lower mean 
value of the aperture in these regions increases its rel- 
ative fluctuations and there are more potential pinning 
sites. As already reported for water in Sec. IIII Bl a buck- 
ling instability occurs and fibers are deformed as they 
get pinned (Fig. [51). Like for water, these deformations 
are localized near the tip of the fiber and move towards 
it as they develop; their amplitude and curvature may 
however become significantly larger than for water be- 
fore depinning takes place (see Fig. [5l for comparison). 
In Fig. [21 for instance, an overhang appears as bump 1 is 
overtaken by 2 as shown by the two upper sets of points in 
the inset. Also, the variation of the radius of curvature 
with time is more complex and not always monotonic 
(Fig. 121) , reflecting the interaction between the different 
loops. The radii of curvature of loops 1 and 2 both be- 
come as low as ~ 2 mm or less: this is lower than the 
radius of loop and, also, much below the radii observed 
with water. 

These larger deformations may result both from the 
larger hydrodynamical forces in the polymer flow due to 
its higher mean viscosity and from the higher depinning 
energies in the narrow part of the wedge. The increased 
value of the hydrodynamical forces is also reflected in the 
faster variations of the distances Sx (see insets of Fig. [5l 
and Fig. [91). 



B. Fiber transport by shear thinning fluids in 
fracture with rough parallel walls. 

1. Aperture field and anisotropy of fracture F3. 

The results of the previous section demonstrate that, 
using a polymer solution, it is possible to transport fibers 
along the full length of model F2, even in the lower 
mean aperture regions. The "wedge like" geometry of 
this model has allowed us to estimate the influence of 
the mean aperture a by varying the distance y at the in- 
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2. Fiber transport in model F3. 



0.2 0.25 
time [s] 

FIG. 9: Variation with time of the radius of the fiber at the 3 
maximum curvature points designated by arrows in Fig. HJd). 
Inset; time variation of the vertical distances 5x between these 
3 points and the fiber tip. Symbols for point 0: (o,»); for point 
1: (□,■) and for point 2: (A, A). 



let of the model where the fiber is injected: however, this 
geometry may have unwanted effects like the sideways 
motions displayed in Fig. 2] and it may mask the influ- 
ence of local aperture fluctuations. The present section 
deals therefore with experiments realized with another 
model (F3) using model rough walls with a similar ge- 
ometry as F2 but with their mean planes parallel (see 
Figs. [TO] for the corresponding aperture map and trans- 
verse mean aperture profile). As mentioned above, the 
fracture walls are two perfectly matching self-affine sur- 
faces which are first put in contact and then pulled apart 
before introducing a small relative shear displacement. 

As shown by several authors 0, H^l , this geometry 
results in an anisotropy of the aperture field and of the 
permeability: flow is easier in the direction normal to 
the shear displacement. In the present experiments, the 
mean flow (x direction) is perpendicular to the shear (y 
direction) and should therefore be oriented in the direc- 
tion of highest permeability. The distribution of the color 
patches in the maps Fig. llOb -b indicates qualitatively 
that the distance over which deviations of the aperture 
from its mean value persist is larger along x than along 
y. For instance, a long continuous vertical string of blue 
patches (corresponding to a high local aperture) is visible 
and is associated with a maximum in the transverse pro- 
file of Fig. [TOb . This suggests a channelization effect re- 
lated to the anisotrop y of the permeability. Recent work 
has indeed shown [23, f26| that, in this configuration, the 
flow field could be described approximately as a set of 
channels of constant hydraulic aperture and parallel to 
the direction x: this allowed, for instance, to predict the 
overall geometry of the displacement fronts between two 
miscible fluids. We shall analyze now how this feature 
influences the transport of flbers in model _F3. 
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FIG. 10: (a) Snapshots of 4 fibers injected into a polymer 
fiow in fracture _F3 at different distances overlaid onto a map 
of the aperture field (mean fiuid velocity U — 215mm/s). 
(b) Snapshot of a fiber coiled up after getting pinned inside 
the fracture {U — 242mm/s). A time sequence of the coiling 
process is displayed in Fig. 1111 The color code corresponds to 
the aperture (red: a < 0.5mm, blue: a > 1.25mm). Field of 
view: 90 x 150 mm. Distances in the vertical x direction are 
counted from the top of the model (see Fig. [TJ. (c): variations 
of the aperture averaged along the distance x parallel to the 
fiow as a function of the transverse distance y. 



In FigurelTOk. snapshots of four fibers injected at differ- 
ent transverse distances y are overlaid onto the aperture 
field of model F3. Three of them are weakly deformed 
and kept moving afterwards down to the lower end of the 
model. One notices that their location does not coincide 
exactly with the zone of highest aperture of the model al- 
though it is not far from it. The fourth fiber (at the left 
on the image) got pinned as it moved out of a high aper- 
ture region into a less open one; the shape of this fiber is 
also more strongly distorted, possibly reflecting a larger 
velocity gradient in the region where the flber moves. In 
model F3, therefore, flber transport is not solely deter- 
mined by the distribution of the aperture even though 
it has a strong influence: other parameters like the local 
gradients of the aperture and slope of the surface or, yet, 
inertia effects may also be of importance. Some of the 
observations on fracture F2 previously suggested similar 
results. 
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3. Irreversible pinning and fiber entanglement in model F3. 
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FIG. 11: Evolution of the fiber as function of time during the 
build-up of a coil for the flow of a polymer solution in fracture 
F3 (see Fig. llOb ). ^From left to right, the time elapsed after 
pinning is respectively 0.63, 1.13, 1.3 and 1.8 s. 

Up to now, we have discussed mostly pinning events 
of finite duration like in Figs. O and [Sj Irreversible pin- 
ning has however also been observed, for instance in the 
narrow parts of fracture F3; this often leads to an entan- 
glement of the fiber around the pinning site and to the 
build up of intertwined loops like in Figs. [Tl] and [TOb . 
Let us discuss in more detail the process leading to the 
configuration displayed in Fig. llOb following pinning in 
a narrower zone (yellow shade region near the top of the 
coil). At first, the fiber gets deformed like in the case of 
transient pinning (see Sections [mB] and lTOI and Figs. [5] 
andlHl). Then, while the tip remains pinned, bump 1 over- 
takes it so that a loop appears at the right of the pinning 
site and move downwards (Fig. [TTk). At the same time 
other parts of the fiber located farther upstream move 
leftwards and build up a second, closed loop (arrow in 
Fig. [TT|) . A key point is the fact that the friction forces 
with the walls in the parts of the fiber located close to the 
pinning site are large enough to prevent hydrodynamic 
forces on fast moving parts of the fiber from dragging 
downward the rest of its length. Tension forces on the 
upstream parts are then quite low so that loops and de- 
formations can easily appear. Finally, both loops move 
down: they get locked on the pinning site and dangle on 
both sides (Fig. [Tib). This further reduces the tension 
on the upstream fiber sections so that new loops start to 
build up (Fig.[nji)- 

V. CONCLUSION 

In the present work, several clear-cut results have been 
obtained experimentally on the conditions for achieving 
the transport of flexible fibers in fractures. First, while 
in smooth fractures the friction of the fibers with the 



walls is negligible and does not influence their motion, 
its influence has been verified to be much stronger for 
rough walls: this may result in pinning and blockage or, 
at least, in a jerky progression instead of a continuous 
one. 

Next, the aperture of the fracture (and particularly its 
ratio to the fiber diameter) is also a key factor, but with 
a more subtle influence than expected. On the one hand, 
in the wedge-shaped rough model F2 with a global trans- 
verse gradient of the mean aperture, fiber transport is, as 
expected, easier in the more open parts of the wedge than 
in narrower ones (however pinning may occur in open re- 
gions). In the rough fracture F3 with globally parallel 
walls and a channelized flow field, fiber motion is also 
easier in high permeability regions. On the other hand, 
the trajectory of the fibers does not coincide fully with 
the paths of highest aperture while pinning is also not 
always localized at minimal aperture points. Other fac- 
tors such as the local slope of the walls and inertia forces 
pushing the fibers towards the solid surfaces must there- 
fore play an important part. In the present experiments, 
the mobility of the fibers appears to increase with the 
mean flow velocity: further experiments will however be 
necessary to determine whether, at high velocities, some 
fibers may be forced into narrow zones by inertia effets 
and get blocked there. 

In order to evaluate the transport of fibers in natu- 
ral or industrial fracture systems of practical interest, 
different configurations of the aperture field need to be 
studied such as, for instance, a channelization perpen- 
dicular to the mean flow. Using our experimental set 
up, this may be investigated by changing the direction of 
the shear displacement between the two complementary 
rough walls with respect to the flow (more precisely by 
comparing for two identical complementary rough walls 
results obtained when the shear displacement u is parallel 
and perpendicular to the mean velocity U). The value of 
Hurst's exponent H for the fracture walls may also influ- 
ence fiber transport by determining the relative height of 
rugosities of different sizes in the plane of the fracture: it 
will be informative to compare the present measurements 
corresponding to H = 0.8 to others achieved for H = 0.5 
(a value encountered in such materials as sandstone and 
sintered glass). 

A third key observation is the enhancement of the mo- 
bility of the fibers when water is replaced by a shear 
thinning fluid at a same flow velocity. For instance, fiber 
propagation along the full length of rough fractures with 
a mean aperture of less than twice the fiber diameter be- 
comes possible. Fiber motion is observed (although in 
a "stop and go" manner) in zones where they got stuck 
in water flows while, in regions of high aperture the dis- 
placement is faster and more continuous. Also, the influ- 
ence of the interaction of the fibers with the walls is much 
reduced and the minimum value of the ratio of the mean 
aperture to the diameter of the fiber for observing trans- 
port is lower. This increased mobility in shear thinning 
flows may result from larger hydrodynamic forces (the 
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effective viscosity of the polymer solutions is generally 
higher than that of water). Other possible explanations 
must however be investigated : for instance, the fiber is 
more frequently localized in the center of the gap between 
the walls than in Newtonian fluids. Also, at high veloci- 
ties, a layer of low viscosity Newtonian fluid appears near 
the walls and may influence the transport of the fibers. 
In order to understand this effect, it will be also impor- 
tant to compare the present results with those obtained 
with Newtonian fluids of different viscosities and other 
non Newtonian fluids such as shear-thinning fluids with 
different characteristic exponents and yield stress or vis- 
coelastic fluids. 

The transport of fibers is also strongly influenced by 
their deformation under the effect of hydrodynamic forces 
which is particularly important during pinning events. In 
the case of intermittent pinning marked by a jerky mo- 
tion, the amplitude of the deformation of the fiber in- 
creases with time and the elastic energy accumulates un- 
til the front tip is released and the fiber recovers roughly 
its initial shape. The amplitude of the deformations and 
the curvatures observed were larger for experiments with 
the polymer solution than for water. In the case of per- 
manent pinning, the amplitude of the deformations may 
increase further: loops appear and pile up on the pinning 
site, finally leading to entanglement. These results show 
qualitatively that, even for very flexible fibers, the influ- 
ence of the elastic modulus on their transport by the flow 
cannot be neglected. 

In future work, systematic investigations will be 
needed to characterize more quantitatively the depen- 



dence of the pinning and entanglement processes on the 
control parameters of the experiments. A first important 
question is the dependence of the motion of the fibers 
on their flexibility and on the value of their diameters 
compared both to the mean fracture aperture a and to 
its standard deviation ctq. Pinning effects also depend 
clearly on the interaction of the fibers and the walls: this 
involves friction between them as well as well as blockage 
of the fiber motion by the wall roughness: their rela- 
tive magnitude will need to be investigated. Overall, as 
shown for instance by Fig. [7| 

Finally, while the present experiments greatly help 
to determine the conditions and characteristics for fiber 
transport in single fractures, further work will be needed 
to understand their transfer from one fracture to another 
in a fracture network. 
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